
function xOut=atanh(xIn)
  x = xIn.v;
  xOut.v=atanh(x);
  d = 1./(1-x.^2);
  xOut.d = elelmult_splastindx(d,xIn.d);

  xOut=class(xOut,'deriv1s');
